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Abstract 

Motivated by the weak-strong coupling expansion |||], we calculate the spec- 
trum of hadrons using a systematic 1/d {d - dimensionality of spacetime) in 



If^ • addition to a strong coupling expansion in (3. The 1/d expansion is pushed to 

the next to leading order in {1/d) for mesons and next to next to leading order 



for baryons. We do the calculation using Wilson fermions with arbitrary r and 
show that doublers decouple from the spectrum only when r is close to the 
Wilson's value r = 1. For these r the spectrum is much closer to the lattice 
results and the phenomenological values than those obtained by using either 
the (nonsystematic) "randomwalk" approximation or the hopping parameter 
expansion. In particular, the value of the nucleon to p - meson mass ratio is 
lowered to ^2arecosh2^ ^i^/'^) ^ 1.48. The result holds even for /? as large as 
5, where the weak-strong coupling expansion is applicable and therefore these 
results are expected to be reasonable. 
It is commonly believed that for low energy physical quantities in QCD, such as hadron 
masses, there is no small expansion parameter, and the theory is "strongly coupled". The 
theory only has an asymptotic weak coupling expansion successfully describes high energy 
quantities but breaks down at low energies. Recently one of us proposed a scheme com- 
bining weak and strong coupling expansion in a "double expansion scheme" [|^. The high 
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energy modes are integrated out using expansion in weak coupling as, the resultant effective 
Lagrangian is then expanded in derivatives and solved using strong coupling expansion in 

A priori it would seem that the domain of applicability for the expansion in j3 will 
have very little chance to intersect with that of the expansion in a. Indeed, if the gauge 
coupling constant g is small, 1/g is large and vice versa. Fortunately, however, a small a{g), 
does not necessarily imply that (3{g) will be large and vice versa. In fact, this scheme of a 
simultaneous weak and strong expansion has been tested on certain solvable low dimensional 
asymptotically free models like Ising chain and d = 2 Gross - Neveu models [0] and the results 
agree very well with the exact values. 

As for QCD, looking at strong coupling expansion results one notices that although 
the asymptotic scaling region is out of reach, the effective weak coupling a near the strong 
coupling radius of convergence is quite small. In the SU{3) YM theory the radius of conver- 
gence (roughening phase transition), is larger than Pmax = ^1 9min ~ ^ [0] which corresponds 
(using naive perturbative RG) to a relatively small effective weak coupling expansion pa- 
rameter aiat = gfat/^'^ ~ O-l- Even taking into account the fact that aiat is a lattice one (not 
the MS, ajj^ ~ 0.3), this corresponds to the values at which perturbation theory is supposed 
to work for energies above 1.5 — 2 GeV. Recently, this fact has been fully understood for 
lattice weak coupling perturbation theory 0. Therefore, there exists a (albeit smaller then 
the one for the Gross - Neveu model) window in the coupling in which both a and (3 are 
small enough to produce a reasonable series. The well known "loop factors" l/(47r)^ in the 
weak coupling expansion parameter a are partly responsible for this, although this generally 
is not sufficient since symmetry group factors tend to reduce the window. Note also that the 
leading order term in a coincides with the conventional "phenomenological" strong coupling 
model, in which the inverse lattice spacing M is limited to values inside the weak-strong 
applicability window. This would seem to be a strong indication that the simultationeous 
weak-strong expansion will be apphcable to QCD. 

The complexity of such a calculation depends on the quantity and the precision one 
would like to achieve. Usually strong coupling series are relatively easy to evaluate to 
very high orders. Consequently the scale M can be chosen in such a way that Pm is just 
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below the strong coupling expansion radius of convergence. Simultaneously, the scale M 
should be sufficiently high or alternatively the relevant energy scale sufficiently low so that 
just a few orders in the derivative expansion are needed to achieve the desired accuracy. 
Consequently, this method is limited to low energy quantities only. Therefore, inside this 
weak-strong coupling window the usual lattice action with renormalized coupling can serve 
as a reasonable effective action 0. The expected accuracy is rather low: up to corrections 
of order ~ 0.3 due to weak coupling expansion. On the top of this we moreover will 
then have to perform the strong coupling expansion. Fortunately it is well known that in 
the hadronic sector (unlike the glueball or pure glue sector in which the above estimates of 
the radius of convergence of the strong coupling were taken) the situation is much better. 
The strong coupling limit (3 = already produces a reasonable spectrum of hadrons. In 
addition, it is known that the next to leading order terms in (3 for the hadronic spectrum 
are numerically very small |10,9[. Even for the window value of /3 = 5 the corrections do 
not exceed 15%. 

However even in this limit QCD is nontrivial and an additional expansion parameter 
should be utilized. This may be the hopping parameter, 1/Nc or 1/d expansions. The 
hopping expansion parameter k ^ 1/2 is defined as a bare coupling in units of lattice spacing 
and therefore cannot be easily related to a physical quantity The 1/Nc = 1/3 has been 
extensively used, but is notoriously difficult to perform beyond the leading order. Indeed, it is 
mostly the nonsystematic random walk approximation P,pi]|,P| that has been used to estimate 
the spectrum at strong coupling. In this approximation the hopping parameter expansion is 
partially summed up, so that quark- lines form "collapsed paths" P,p!0|. Although it seems to 
be superior than the simple hopping parameter expansion, one does not find a controllable 
expansion parameter within this approximation. Moreover, the results of this approximation 
as well as those of the hopping parameter expansion were rather discouraging. Although 
the ordering of the lowest hadronic states is correct, some mass ratios are grotesque. An 
especially bad example is the nucleon to p - meson mass ratio which is about 2.2 instead of 
the phenomenological 1.2 or (at [3 = 5) lattice MC simulation value of ~ 1.4. 

In this paper we use the 1/d = 1/4 expansion to calculate the hadron mass spectrum, 
which allows managable higher order calculations. This was ffist applied to Yang-Mills 
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theories for staggered fermions in The results for the nucleon to p mass ratio is ~ 1.7, 



which, although better the previously mentioned strong coupling result is still so different 
from the phenomenological value that it cannot be accounted for by next to leading order 
corrections. Moreover, the interpretation of particle spectrum for staggered fermions is by 
no means straightforward. We, therefore, shall consider Wilson fermions which allows us to 
discuss splitting due to three flavors and for which the interpretation is straightforward. 

In view of these above results for the spectrum, one of the following should be valid: 
(a) something is nevertheless wrong with the argument for the existence of the weak-strong 
coupling window and the quenched lattice results are not precise enough. The spectrum at 
/? = 5 is indeed very different from the observed experimental one because the continuum 
limit is still far from this point or (b) the random walk and the hopping expansion (and to 
a smaller degree the staggered fermion) results are inaccurate. We show in this paper that 
when a systematic 1/d expansion is applied to the Wilson action within the weak-strong 
window, a spectrum is obtained which is in agreement with the above lattice MC results 
within the expected accuracy of the expansion. 

We now fix notation and outline the formalism, which is well-known, focusing on the 
differences with random walk approach. No details of higher order calculations will be 
given. The standard lattice Wilson action is 

c plaquettes x,fi x 

where J^^(x) = ^/'^(a; + yu)P+?/'f (x) and = (r ± 7^)/2. f/'^^ is the usual compact gauge 
field on the lattice and the script letters run over A^^ colors. The lower latin letters runs over 
Nf flavors. First, we shall limit ourselves to the (3 = limit and then discuss the effects of 
finite /?. 

Integration over the gauge fields U give in the leading order H], 



det J;f^(x) + det J-^^(x) 



+ ... (2) 



where the determinant is over the color indices. The . . . indicates a finite number of terms 
which contribute to higher order terms in 1/d. Considering first the mesonic sector, we 
note that although we have introduced different flavors, for the leading order in 1/d, they 



will not play a role and we shall suppress them for now. We then introduce the mesonic 
fields through M'^(x) = -^sM ip'^T'^tfj-^. The "channel" index A runs from to 15 and 



normalization of matrices F is chosen so that Tr r^^F^ = 6^^ ||T0| . In terms of these fields 
the mesonic part of the action becomes 

xyfj. 

where unless otherwised stated summation over repeated channel indices will be understood. 
We now integrate over the fermion fields. This is done by first introducing auxiliary fields 



conjugate to M [|12 



omes 



I VM^e^v\NcY.M\x)M\x)-^Y.J^\^)DlB{x-y)M^{y)\ (4) 
Then the remaining gaussian integral over fermionic fields can be done: 
Z, = j exp ^ M''{x)DfB{x - y)M''{y) - N,J2 Tr^ log (T^M'^ix) + 2m) | 

= J PA<^e--^[^(")] (5) 



where fn = m/ \f2d. The factor ^/d was introduced in the mass to facilitate the l/ci expansion 
||1 1|| . The functional ^[A^(a;)] of hadronic fields is an effective hadronic action describing 
dynamics of the color invariant "basic" fields only. In the meson sector, which is being 
considered now, these fields interpolate between the pseudoscalar and vector mesons. They 
correspond to the lowest energy states of the naive quark model. Later on baryonic fields 
interpolating between the octet and decouplet [N and A) fields will be introduced. 
The quadratic part of the mesonic action to lowest order is then 

= -^Y.^\^)GtB{^-v)M\y) (6) 

xy 

where 

G~ab{^ - y) = Dfj^{x - y) \ ^ (7) 



and Ao = + \/rn + 1 — comes from the solution of the gap equation [11|. To find 
the mass spectrum of the theory, we need to find the zeros of in momentum space. 
This would seem to be difficult, but note that G^^ and k = AqG"^/^ have the same zero 



eigenvalues, as long as D does not vanish and finding the zeros of the latter quantity is fairly 
straightforward. 

There are two coupled channels for the mesonic sector. First, the pseudoscalar couples 
with the time component of the axial vector giving the mass term for the pion. In this 
channel. 



/ An - (1 + r^) + (1 - coshm) -i^ sinhm 



G-'D\ 



—i^ sinhm Aq + (1 — r^) — — ^^^J -* coshm 



(8) 



where we included the correct powers of l/ci and have taken the d- momentum to be 
(0, ...,0,im). Within the 1/d expansion, we can generically write 

cosh[m] = xci + ?/ + 0(1). (9) 

where, of course, y cannot be determined as yet since higher orders in l/rf terms have not 
been included. Using eq.(H) and expanding the eigenvalue equation of the matrix eq.(7) in 
orders of 1/d, we find two solutions. The one which is finite for all relevant r, is 

X. = {{XI - r'){l + r') - [(1 - r^? + ^r\Xl - r^] '^'] (10) 

and determines the mass of the mass of the pion. The second eigenvalue Xd is nonzero and 
describes a doubler (it is a bound state of two fermionic doublers at the opposite corners of 
the Brillouin zone): 

2rf(l + r2) 

coshmrf = — — (11) 

(1 — r^)^ 

One then requires that the pion be massless for all r, and sets = and = 1. This, 
to lowest order, determines the trajectory in parameter space relating bare mass to r: Aq = 
1 + r^. 

In the p meson channel the corresponding matrix (on the trajectory) is: 

Aq - (1 + r^) + i (3 + - (1 + r^) coshm) -i^ sinhm 

sinhm Ao(l — r^) + ^[—3 + — (1 + r^) coshm] 

(12) 
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(13) 



where we have considered only one spatial component of the vector field interpolating p 
mesons. On this trajectory we write cosh mp = Xpd + yp + 0{l/d) and obtain 

— (l + r^)xp —2irXp 
2irXp 2 — (1 + r'^)xp 
to the leading order in 1/d. There are two eigenvalues to this matrix. The first is Xp = 0. 
This means that the p meson's mass does not have a "natural" order of log d, but is , in fact, 
smaller - just of order 1 . It is, however, inconsistent to to determine i/p by solving equation 
for vanishing determinant of eq.(|12|). The 1/d corrections to this matrix, considered in the 
following, must be taken into account. Other channels, scalar, tensor, etc, contain doublers 
only. 

Genericly, in any next to leading order calculation in 1/d there are two types of contri- 
butions The first is the "tree" contribution which arises from additional terms in the 



integral over gauge fields eq.(^, while the second is the one loop diagrams involving the 
propagator and vertices of the leading effective action eq.(^). Keeping the pion mass zero, 
the "next to leading" contribution to the p meson's mass yp (which is actually the leading 
since Xp = 0) is: 

This value is consistent with the r = result obtained in for staggered fermions. Note 



that for r = 1 the p mass is significantly larger then in the random walk approximation. 

The purpose of introducing the chiral symmetry breaking mass and Wilson's terms was 
to remove doublers. The value of r should, in principle, be optimized in such a way that on 
the one hand doublers do not interfere with physical particles and, on the other hand, the 
chiral symmetry is minimally violated. As is well known, setting the pion mass to zero does 
not mean that the chiral symmetry is somehow restored on the trajectory. It just means 
that we are situated on the spontaneous parity breaking phase transition line. On Fig.l we 
show the r dependence of various doublers masses compared to the p meson mass. We can 
see that the doublers are significantly heavier then p mesons only near Wilson's value of 
r = 1, where they all become infinitely heavy. Therefore we conclude that there is no great 



advantage to work with r < 1 contrary to some lattice [jT3| and random walk approximation 
fSll results in which doublers were not considered. 



Turning our attention to the baryonic sector, we note that the general expression for the 
exponent of the baryon masses contains, half integer as well as integer powers of 1/d: 



S'^ + x'd + yS'^ + y' + 0{l/^) (15) 



We calculated x, x', y and y' . To the leading order the A mass on the = trajectory is 



2v/ 2(l + r ^^3/2 

:i+r) 



XA = TTT-Ts (16) 



Notice that for r = 0, xa = 2\/2 which reproduces the result obtained in |TT[| for staggered 
fermions calculation. The nucleon is degenerate with A up to the order of 1/d we have 
considered when r = 1. The formulas for the higher order corrections to the baryon masses 
for other values of r are cumbersome and will be given elswhere. Instead, for r = 1 

1 1 , 29 

^A = -4, 2/A = --, yA = Y^ (17) 

We also calculated the leading order /3 corrections to the p and baryons. These types of 



corrections has been studied for staggered fermions in ||Tl|] and within the random walk 
approximation scheme in . They vanish for r = 1 and are small for other values of r near 
the window range. 

To summarize, we have systematically studied the spectrum of the Wilson action with 
arbitrary r using strong coupling and the 1/d expansions. The Lagrangian is considered as a 
phenomenological low energy effective Lagrangian for the values of (3 at which the expansion 
is still reasonable. It turns out that the doublers decouple from the physical spectrum only 
when r is quite close to the Wilson's value r = 1. The results of the systematic l/d expansion 
for the nucleon to p mass ratio is ^2^^008112^ 0{l/d) ~ 1.48. This value is within the range 
of the next correction of the Monte Carlo results. 

We now compare the 1/d results for mesons with those obtained within the random walk 
approximation in = 4. To understand the difference with the random walk approach, we 
have extended the random walk calculation to arbitrary d and obtained the mass of p 
meson for r = 1 as: 

2d-l 3 3 3 

eo.h™, = — = 2-- + ^-^ + ... (18) 
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The first term in this expansion coincides with the systematic 1/d expansion. There is 
however no reason to expect that the next order term in 1/d will resemble that in eq.(^). 
A priori this is not obvious since at the special value r = 1, many contributions to the next 
to leading order term vanish. This is due to the well-known result from the random walk 
approximation that to this order all the propagators travel along a single direction. As 



such, we will necessarily encounter product of two projectors P^P/^ = 0. However if we go 
to the next to next to leading order there will certainly be many contributions which will 
not vanish for the simple reason that certain diagrams will not "linear", but rather "planar". 

The coefficient in front of l/rf in eq.(0) is negative and large numerically for 4 = 4. This 
leads to significant underestimate for the p mass and consequently for the overestimation 
of rriN/mp. Let us emphasize that this term cannot be taken seriously at this point since 
corrections to the mass at this order in l/rf have not been done as yet. Indeed, it would be 
very interesting to calculate this next {1/d for coshmp) order term to obtain better estimate 
of the p mass. 

This feature of random walk approximation does not carry over to baryon sector, however. 
The corresponding expansion for the nucleon mass in powers of 1/d is: 



The first two terms of these now coincide with our systematic 1/d expansion results. When 
a similar expansion is done for the Delta mass, we find that for some unknown reason all 
four terms now coincide with eq.(0). Once again, however, any terms which is of higher 
order than d in eq. ([19|) are unreliable, since they will be almost certainly be changed by 
higher order 1/d corrections. In particular, we see that small splitting between the nucleon 
and A found in random walk approach is due these unreliable higher order terms. Within 
the systematic 1/d expansion, the Delta and Nucleon are degenerate. 

Of course there are numerous other corrections to the weak-strong approximation scheme. 
Here we discuss few of the many which have not been calculated. Our previous discussion 
revealed the peculiar fact most of the contributions to the next to leading order terms in 1/d 
vanish for r = 1. We can ask the following question: What are the corrections that do not 
vanish? In particular, we note that to the next to leading order there is no splitting between 
Goldstone bosons and fiavour singlet - the r] particles. As is well known, this splitting is 



due to an anomaly. The plaquette correction also does not lead to the splitting. This is 
easy to understand. The arguments of Frohlich and King are applicable to the 1/d 



expansion. We therefore expect these anomaly effects to appear only at very high orders in 
1/d or (3, when a diagram which spreads in all four directions can be constructed. Another 
possibility is that the main mechanism for the splitting is not due to these corrections but 
rather to direct anomaly breaking terms which are proportional to a^- These appear due to 
the presence of instantons at energies higher then the scale M. 

The actual splitting between A and the nucleon is also probably due to higher order terms 
in the effective Lagrangian. Note that in lattice simulations at relatively small /? (/5 ~ 5) the 
splitting is also smaller then the phenomenological values. The presence of higher dimension 
terms are also crucial for two other purposes: restoration of the Lorentz invariance and the 
chiral symmetry. Chiral symmetry is explicitly broken by the mass and Wilson terms. As we 
have mentioned, the vanishing of the pion mass is not sufficient tp restore chiral symmetry. 
Instead it is simply a signal criticality with respect to a discrete symmetry. It would be 
therefore be interesting to investigate whether the chiral properties are gradually restored 
once higher dimensional operators are introduced. For example the pion scattering at small 
momenta is nonvanishing without higher dimensional terms. It is reasonable to expect 
that with the inclusion of the term due to next order correction in of the derivitive 

expansion of the effective action the correct zero momentum limit at least will be recovered. 
Existing results, which are quite scarce within the random walk approximation scheme, in 
which only part of the higher dimensional (improvement) operators are considered did not 
address this question. 
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Fig.l Masses of nucleon, p meson, mesonic and baryonic doublers as function of Wilson 
parameter r. 
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